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EFFECT  OF  NONUNIFORM  SIZE  ON  INTERNAL 
STRESSES  IN  A  RAPID,  SIMPLE  SHEAR  FLOW  OF 
GRANULAR  MATERIALS 
Part  1.  Two  Grain  Sizes 

Hay  ley  H.  Shen 

INTRODUCTION 

Many  industrial  and  geophysical  processes  involve  the  transport  of  solids  in  granular  form.  Im¬ 
portant  example.,  include  flows  in  hoppers,  chutes  and  slurry  pipelines,  sediment  transport  in  rivers, 
snow  avalanches  and  debris  flows.  All  of  these  examples  deal  with  moving  particles  mixed  with  air 
or  a  liquid,  and  contain  solids  that  can  range  in  size  from  boulders  to  a  very  fine  dust.  At  present 
all  quantitative  descriptions  of  such  flows  are  largely  empirical  or  semi-empirical  (Durand  1953, 

Newitt  et  al.  1955,  Voellmy  1955,  Kennedy  and  Brooks  1965,  Bosley  et  al.  1969,  Roberts  1969). 

Recent  theories  describing  the  constitutive  relationships  for  rapidly  sheared  granular  mixtures 
have  revealed  many  insights  into  the  properties  of  the  flows  of  such  mixtures  (Jenkins  and  Satake 
1983).  These  theories  have  yielded  stress-strain  rate  relationships,  and  they  consider  binary  colli¬ 
sions  between  the  granular  particles  to  be  the  major  stress-generating  mechanism. 

Comparisons  between  these  theoretical  results  and  experimental  data  obtained  with  various 
materials  (Bagnold  1954,  Savage  and  Sayed  1980,  Savage  and  Sayed  1984)  were  made  by  Shen 
and  Ackermann  (1982).  The  experimental  data  were  obtained  by  shearing  dry  or  neutrally  buoy¬ 
ant  spherical  particles  of  uniform  diameter  in  a  couette  device.  The  comparisons  show  that  current 
theoretical  developments  are  capable  of  not  only  revealing  qualitative  insights  into  granular  flow, 
but  are  also  capable  of  predicting  quantitatively  measured  results  quite  accurately. 

The  major  limitation  to  applying  the  theories  to  real  granular  flows  is  that  the  constitutive  re¬ 
lationships  have  been  developed  for  granular  mixtures  of  identical  size.  Almost  all  real  granular 
flows,  however,  consist  of  a  wide  range  of  particle  sizes.  The  best  an  engineer  can  now  do  in  anal¬ 
yzing  a  granular  flow  is  to  use  a  somewhat  arbitrary  “characteristic”  size  to  establish  empirical 
formulas. 

Experiments  have  shown  (Gilbert  1914,  Durand  1953,  Savage  and  Sayed  1984)  that  grain 
size  distribution  has  a  significant  effect  on  flow  properties.  Consider,  for  example,  the  data 
obtained  in  an  experiment  in  which  two  mixtures  were  identical  in  all  respects  except  one;  the 
first  mixture  was  formed  from  solids  having  a  grain  diameter  of  1 .32  mm,  while  the  other  mixture 
was  formed  from  two  particle  sizes  (0.6  mm  and  1 .65  mm)  combined  to  produce  a  mean  particle 
diameter  of  1.34  mm,  nearly  equal  to  that  of  the  first  mixture  (Savage  and  Sayed  1984).  The 
mixture  with  a  uniform  particle  size  generated  stresses  that  were  almost  five  times  higher  than 
those  generated  with  a  mixture  of  two  sizes,  even  though  the  mean  diameters  were  similar. 

In  Part  I  of  this  study  the  constitutive  equations  are  developed  to  describe  the  stress-strain  rate 
relationships  for  rapidly  sheared  mixtures  in  which  the  granular  solids  have  two  different  particle 
diameters.  The  extension  to  a  complete  distribution  of  particle  sizes  will  appear  in  a  future  paper. 

In  Pari  2  of  this  study,  similar  results  are  obtained  for  multiple  grain  size  mixtures  (Shen  1 985). 

As  was  done  by  Jenkins  and  Savage  ( 1983)  and  Shen  and  Ackermann  (1982),  the  binary  collisions 
between  particles  are  considered  the  major  internal  momentum  transfer  mechanism  that  gives  rise 
to  stresses.  Comparisons  between  the  theoretical  results  and  experimental  data  show  good  agree¬ 
ment. 
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Figure  1.  Shear  flow  of  a  mixture  of  spheres.  Figure  2.  A  control  volume  in  a  granu¬ 

lar  flow  of  uniform  spheres. 

MODEL  OF  A  RAPIDLY  SHEARED  GRANULAR  FLOW 


Consider  a  simple  shear  flow  of  a  granular  mixture  that  is  steady  and  uniform  as  shown  in  Fig¬ 
ure  1 .  Assume  that  the  solid  particles  are  all  made  of  the  same  material  but  have  two  sizes.  Let 
Dl  and  Ds  represent  the  diameters  of  the  large  and  small  particles,  respectively,  and  ps  equal  the 
density  of  these  particles.  Let  e  be  the  restitution  coefficient  and  p  the  kinetic  friction  coefficient 
of  the  solids.  The  interstitial  fluid  effect  will  be  neglected  at  the  moment.  The  shear  flow  has  a 
linear  velocity  profile  with  u(x2 )  representing  the  mean  local  velocity  at  coordinate  location  x2 . 

It  is  also  assumed  that  the  body  forces  are  negligible.  When  the  velocity  gradient  is  constant  and 
the  body  forces  negligible,  the  concentrations  of  both  types  of  particles  can  be  assumed  to  be 
uniform  throughout  the  flow  field. 

When  particles  are  sufficiently  numerous  and  the  shearing  rate  is  sufficiently  high,  the  differ¬ 
ences  in  the  mean  motion  of  adjacent  particles  will  trigger  collisions  among  the  surrounding  parti¬ 
cles.  These  collisions  will  in  turn  produce  additional  collisions  and  random  fluctuations  of  all  other 
particles  through  the  flow  field.  A  condition  is  therefore  postulated  in  which  every  particle  has  a 
mean  velocity  u(x2 )  and  a  continually  fluctuating  velocity  component  superimposed  on  the  mean 
motion.  This  type  of  flow  field  is  consistent  with  flow  conditions  of  a  granular  mixture  postulated 
in  previous  analyses  (Kanatani  1979,  Ogawa  et  al.  1980,  Shen  and  Ackermann  1982,  Jenkins  and 
Savage  1983)  and  has  been  observed  in  experimental  studies.*  The  collisions  create  a  momentum 
transfer  between  particles.  The  rate  of  the  net  momentum  transfer  across  any  unit  surface  area 
can  then  be  interpreted  as  creating  both  shear  and  normal  stresses  on  that  area. 

As  the  frequency  of  collisions  rises,  the  resulting  momentum  transfer  eventually  dominates  any 
other  mechanism  that  might  also  generate  stresses.  When  stresses  are  generated  predominantly  by 
momentum  transfer  from  collisions,  the  granular  flow  is  described  as  being  in  the  “grain  inertia" 
region  (Bagnold  1954).  As  in  previous  analyses  (Kanatani  1979,  Ogawa  et  al.  1980,  Shen  and 
Ackermann  1982,  Jenkins  and  Savage  1983)  the  basic  assumptions  are:  binary  collisions  are  the 
dominating  mechanism  for  stress  generation,  viscous  stresses  generated  by  fluid-fluid  shearing  are 
negligible,  and  the  mass  transfer  and  associated  momentum  transfer  are  negligible  across  any  sur¬ 
faces  parallel  to  the  How  direction.  For  mixtures  of  granular  solids  having  water  or  air  as  the  inter¬ 
stitial  fluid,  the  concentrations  that  will  make  the  above  assumptions  valid  should  be  more  than 
2(Y/r  by  volume  and  the  shearing  rate  more  than  approximately  10  s"1.  These  criteria  have  been 
obtained  from  the  experimental  work  done  by  Bagnold  (1954)  and  Savage  and  Sayed  (1980). 

The  stresses  on  a  unit  surface  are  defined  as  the  total  rate  of  momentum  transfer  across  a  unit 
surface.  Under  the  above  assumptions,  this  momentum  transfer  is  caused  by  the  binary  collisions 
between  particles  located  on  the  unit  surface  and  particles  coming  from  outside  a  unit  volume 


Persona!  communication  with  N.L.  Ackermann,  C'larkson  University,  Potsdam,  N.Y.,  1984. 


bounded  by  that  unit  surface.  For  illustration, 
consider  a  granular  flow  of  single-size  particles. 

Let  p2  be  the  number  of  particles  on  a  unit 
surface  perpendicular  to  the  x2 -direction,  as 
shown  by  the  top  surface  in  Figure  2.  Let  / 
be  the  frequency  of  binary  collisions  between 
a  particle  on  that  surface  and  the  particles  ex¬ 
terior  to  that  surface.  Let  A A/,  be  the  average 
momentum  transferred  in  the  .v,  direction  to 
the  surface  particles  after  each  collision.  Tire 
shear  stress  on  that  surface  t2  i  can  therefore 
be  computed  as 

*21  =  P:  '  1  '  At/)  .  ( 1 ) 

The  normal  stress  on  the  surface  t2 2  is  computed  similarly,  with  AW,  being  replaced  by  AW2 , 
where  AW2  is  the  average  momentum  transferred  in  the  ,v2 -direction.  The  above  equation  was 
first  established  by  Bagnold  ( 1954).  The  idea  used  to  develop  eq  1  can  be  extended  to  flows  with 
particles  that  are  not  uniform,  with  the  right  side  of  eq  1  being  appropriately  generalized  to  incorp¬ 
orate  combinations  of  different  particles. 

FORMULATION  OF  STRESSES  IN  A  TWO-SIZE  MIXTURE 

Consider  a  unit  volume  in  the  shear  flow,  as  shown  in  Figure  3.  On  the  top  surface  there  are 
four  different  kinds  of  collisions:  a  large  exterior  particle  colliding  with  a  large  surface  particle, 
a  large  exterior  particle  colliding  with  a  small  surface  particle,  a  small  exterior  particle  colliding 
with  a  large  surface  particle,  and  a  small  exterior  particle  colliding  with  a  small  surface  particle. 

Any  collision  between  the  surface  and  exterior  particles  results  in  a  momentum  transfer  across 
the  surface.  This  momentum  transfer  is  attributable  to  the  differences  in  the  initial  momentums 
of  the  particles.  Tire  initial  relative  momentum  of  colliding  particles  depends  on  not  only  the  dif¬ 
ference  between  their  mean  velocities,  but  also  the  difference  in  their  fluctuating  velocities. 

The  following  notations  will  be  used  in  the  ensuing  discussion: 

;  =  1  or  2,  indicating  the  v,  -  01  v2 -direction. 

AWLSj  =  average  momentum  AW  transferred  in  the  r  edirection  from  the 
collisions  between  a  large  surface  particle  and  a  small  exterior 
particle. 

AWsl 1  =  average  momentum  transferred  in  the  r  edirection  from  the  colli¬ 
sions  between  a  small  surface  particle  and  a  large  exterior  particle. 

AV/LLj  and  AWssj  =  average  momentum  transferred  in  the  .r  edirection  from  collisions 

between  a  surface  particle  and  an  exterior  particle  of  the  same  kind. 

A'J  s  =  frequency  of  collisions  between  a  single  large  surface  particle  and 
any  of  the  surrounding  small  particles. 

A'^t  =  frequency  of  collisions  between  a  single  small  surface  particle  and 
any  of  the  surrounding  large  particles. 

A 1  j  and  .V£s  =  frequencies  of  collisions  between  a  single  surface  particle  and  the 
surrounding  particles  that  are  of  the  same  kind. 

Pi  -P s  =  number  of  large  or  small  particles  on  a  unit  surface. 
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Figure  3.  A  control  volume  in  a 
granular  flow  of  mixed  spheres. 


If  we  assume  that  the  concentration  and  random  motion  of  both  types  of  particles  are  uniform, 
the  collision  rate  between  a  single  surface  particle  and  the  exterior  particles  should  be  half  of  the 
collision  rate  between  a  surface  particle  and  all  the  surrounding  particles.  Hence,  Vs  V  jV|l/2, 
AlL/2  and  jV|s/2  indicate  the  four  collision  rates  of  the  four  different  types  of  collisions  between 
a  surface  particle  and  the  exterior  particles. 

Using  the  notation  defined  above,  we  can  formulate  the  stresses  as 
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To  determine  the  collision  frequency,  the  energy  balance  equation  will  be  utilized.  First,  a 
concept  of  internal  energy  must  be  established.  The  collisions  that  take  place  in  a  rapidly  sheared 
granular  flow  cause  the  particles  to  move  randomly.  This  random  motion  can  be  represented  by 
velocity  fluctuations  about  the  mean  velocity.  For  the  case  of  uniform  particles,  let  v  be  the 
fluctuation  velocity  about  the  mean  velocity  for  any  particle.  The  internal  energy  per  unit  volume 
contained  in  the  solids  is  then  defined  as  !4  C  psv'2 ,  where  C  is  the  volume  concentration  of  solids. 
In  a  steady,  simple  shear  flow  with  a  uniform  shearing  rate,  the  energy  balance  is  written  as  (Shen 
and  Ackermann  1982) 


L114  i  i  r-i  j  > 

Tu^=N.F.t' 


On  the  left  side  of  eq  4  is  the  rate  of  work  spent  on  increasing  the  internal  energy  of  a  unit  volume 
in  the  flow.  On  the  right  side  of  eq  4  is  the  rate  of  internal  energy  dissipation  in  a  unit  volume  be¬ 
cause  of  collisions.  In  eq  4  N  denotes  the  number  of  particles  in  a  unit  volume,  F  the  collision  fre¬ 
quency  per  particle  and  E  the  average  energy  dissipation  per  particle  per  collision. 

For  a  mixture  of  two  different  sizes,  the  energy  equation  described  by  eq  4  should  be  general¬ 
ized  as  follows: 


r21  Jx  ~  Vs  Vs  +  Vl  ^  LL  +  Vs  Vs 


where  Vs  is  the  rate  of  collisions  between  all  pairs  of  large  and  small  particles  in  a  unit  volume, 
while  Vl  and  Arss  are  the  rates  of  binary  collisions  between  like  particles  in  a  unit  volume,  and 
Vs  is  the  average  energy  loss  in  a  pair  of  colliding  large  and  small  particles,  while  A'Ll  and  /TSs 
are  the  average  energy  losses  in  a  pair  of  colliding  like  particles.  After  substituting  eq  2  into  eq  5, 
we  obtain 


'LL!  2 


'SLl  2 


'ssi  2 


Vl  ^  LL  +  'V.S  Vs  +  Vs  Vss  • 


The  parameters  that  appear  in  the  above  equation  are  functions  of  the  flow  condition  and  material 


properties.  These  functions  can  be  summarized  as  follows  in  eq  7-1 1 : 


Pi  or  S  PLorS^L’^S’^L'^s)  (?) 

where  CL  and  Cs  are  the  concentrations  of  the  large  and  small  particles  and  Dh  and  Ds  are  their 
diameters; 

A*hq  =  Kq  (Cl ,  Cs,  Dl,  Ds,  v'h ,  v's)  ;  P  and  Q  =  L  or  5  (8) 

where  i'j_  and  v's  are  the  fluctuation  velocity  of  large  and  small  particles,  respectively; 

^vq^Nvq(ClXs,Dl,Ds,v'l,V's)  -  PandQ  =  LorS  (9) 

AM  rot  =  MPQi  ,  v'L  ,v's,Dl,Ds,  Cl  ,Cs,e,  h,  fis)  ;  P  and  Q  =  L  or  5,  (10) 

/  =  1  or  2 


and 


/■'hq  =  /•-pq('^-  .  '-’l  .  •'s-  dl ’ ^s*  e-  /h  Ps) ;  Z’  and  Q  =  L  or  5. 


(11) 


In  eq  7-1 1  the  parameters  du/dx2 ,  DL,  Ds,  CL,  Cs,ps,  e  and  p  are  specified.  The  unknown  vari¬ 
ables  are  the  fluctuation  velocities  v'L  and  v's. 

After  the  functional  relationships  described  in  eq  7-1 1  are  established,  the  fluctuation  speeds 
r'L  and  i  ^  can  be  determined  using  eq  6  and  by  assuming  an  equipartition  of  energy.  The  equi- 
partition  of  energy  is  described  as 


(12) 


which  says  that  the  internal  energy  (the  “random  fluctuation  energy”)  of  the  small  and  large 
particles  is  the  same.  As  proved  in  many  texts  on  the  kinetic  theory  of  gases  (Kennard  1938, 
Chapman  and  Cowling  1970)  the  random  fluctuations  of  different  particles  that  form  an 
equilibrium  mixture  contain  the  same  amount  of  fluctuation  energy.  In  a  steady,  constant 
velocity  gradient,  rectilinear  granular  flow,  the  mechanism  that  generates  the  random  fluctua¬ 
tion  energy,  is  uniform  everywhere.  Hence,  a  random  energy  equilibrium  state  can  be  assumed. 

The  following  section  is  devoted  to  explicitly  quantifying  eq  7-1 1 . 


MICROSCOPIC  KINEMATICS  AND  DYNAMICS 

Number  density  of  large  and  small  spheres  on  a  unit  surface -pL,  Ps 

Consider  a  unit  surface  area  perpendicular  to  the  x2  -direction.  This  unit  surface  cuts  through 
numerous  large  and  small  particles  at  various  depths  in  the  particles.  If  we  assume  a  completely 
random  situation,  the  average  area  cut  by  the  unit  surface  can  be  obtained  by 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
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Shen,  Hayley  H. 

Effect  of  nonuniform  size  on  internal  stresses  in  a 
rapid,  simple  shear  flow  of  granular  materials:  Part  1. 

Two  grain  sizes  /  by  Hayley  H.  Shen.  Hanover,  N.H.:  Cold 
Regions  Research  and  Engineering  Laboratory;  Springfield, 

Va . :  available  from  National  Technical  Information  Ser¬ 
vice,  1985. 

iv,  25  p.,  illus.;  28  cm.  (  CRREL  Report  85-2.  ) 

Bibliography:  p.  13. 

1.  Flow.  2.  Granular  flow.  3.  Slurries.  I.  United 
States.  Army.  Corps  of  Engineers.  II.  Cold  Regions  Research 
and  Engineering  Laboratory,  Hanover,  N.H.  III.  Series: 
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Substituting  eq  B1-B7  into  eq  B5  yields 


p  __  nDL  v_l_+vs  rl-e2  ^  a<(  1  +e)  /i2(l+e)2 
as  "Ps  6  D.  .  1  4  n  '  4 

us 


(B8) 


APPENDIX  B:  AVERAGE  ENERGY  LOSS  IN  A  PAIR  OF  COLLIDING  LARGE  AND  SMALL 
SPHERES 


Consider  a  pair  of  large  and  small  spheres  about  to  collide,  as  shown  in  Figure  B1 .  The  total 
energy  loss  because  of  a  collision  is  independent  of  the  inertial  reference  frame;  hence  a  reference 
frame  is  chosen  to  move  with  the  center  of  mass  of  these  two  spheres.  If  we  assume  that  only 
fluctuation  components  just  prior  to  collision  are  important  in  this  computation,  the  velocities  of 
the  large  and  small  spheres  in  the  normal  and  parallel  directions  of  contact  are  related  by  (Green¬ 
wood  1965) 


1'ns 


(Bl) 


Figure  Bl.  A  collision  between  large  and 
small  spheres,  center  of  mass  system. 


The  average  momentums  transferred  from  a  small  sphere  to  a  large  sphere  after  the  collision 


are 


jt/2  3tt/2 


—  2  r  ~  r  n/)? 

^LSI  J  d<t>  J  VNL(N  +  ^P)  ‘  ide 


(A6) 


and 


nil 


n/2 

3n/2 


AWLS2  =  “l/  J  VNL  (A'  +  pf>)  ’fdO  . 

0  7t/2 


(A7) 


The  term  juF  in  eq  A7  will  disappear  under  the  assumption  that  e'L  or  s  »  (du/dx2)[(DL+Ds)/2] . 
This  argument  was  also  adopted  by  Shen  and  Ackermann  (1982).  Substituting  eq  A1-A5  into  eq 
A6  and  substituting  eq  A1-A5  into  eq  A7,  we  obtain 


nDi 


2D\ 


^ls,  =Ps-fL(He)(0.053  +  0.081  p)|^(^_+  s)-^ 


(A8) 


and 


A1-  nDl  4  0+e)Os  ^ - - 

^LS2  =  -Ps  ~r  ^  HTTTvr  +  ‘-r2  • 


-•  6  n2  +  Z)g  v  s  1  't 


(A9) 


In  a  collision  process  only  impulsive  forces  are  involved;  hence  the  momentum  transferred  from  a 
large  sphere  to  a  small  sphere  equals  that  from  a  small  sphere  to  a  large  sphere,  or 


^LSl  or  2  ~  A^SLl 


or  2 


(A10) 
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APPENDIX  A:  AVERAGE  MOMENTUM  TRANSFER  BETWEEN  A  PAIR  OF  COLLIDING 
LARGE  .AND  SMALL  SPHERES 


Consider  a  pair  of  large  and  small  spheres  about  to  collide,  as  shown  in  Figure  A1 .  The  refer¬ 
ence  frame  is  fixed  to  move  with  the  center  of  the  large  sphere  at  a  velocity  equal  to  the  velocity 
of  the  large  sphere  just  prior  to  collision.  For  the  large  sphere  in  this  reference  frame,  the  post¬ 
collision  velocities,  as  indicated  by  *,  in  the  directions  normal  and  parallel  to  the  contact  surface, 
are  (Greenwood  1965) 


(M-e)Ds 
D\  +  I)i 


(l+e)£>g 

^  d3l+d3s  l> 


where  rNS  is  the  normal  component  of  the  small  particle’s  velocity.  An  assumption  is  made  in 
deriving  rfL  such  that  the  impulsive  frictional  force  is  proportional  to  the  impulsive  normal  force, 
with  proportionality  p  being  the  kinetic  frictional  coefficient. 

If  we  follow  the  argument  used  in  deriving  stresses  for  a  mixture  of  uniform-sized  spheres 
(Shen  and  Ackermann  1982)  when  computing  the  average  momentum  transfer  in  the  at, -direction, 
the  velocity  of  the  small  particle  is  assumed  to  be 

v- s  =  <-~i  +s)  ^  cos0i  .  (A2) 

2  dx2 

When  computing  the  average  momentum  transfer  in  the  x2 -direction,  it  is  assumed  that 


?s  =  vi's  +  1’l  (cosa/V  +  sinaP )  ;  U<a<y  .  (A3 

The  two  unit  vectors  IV  and  P  indicate  the  normal  and  parallel  directions  of  collisions.  They  are 
quantified  as 

N  =  -  sin0(cos#  /  -  sin#  k)  -  cos0  /  (A4 


f — FT  .  i,  .  r*  sin0cos0cosf?/  ,  sin2  0cos0  sinOA 

P  =  Vcos2  0+stn2  0sin2  6  [ / - ? - f-* -A+  — — 5-— • 

1  cosJ0+stn20sin20  cos2  0+sin2  0sin2  0 1 


f-'igurc  Al.  A  collision  between  large  and 
small  spheres. 
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be  included  at  a  crude  level  as  was  done  by  Shen  and  Ackermann  (1982),  in  which  the  energy  dis¬ 
sipation  attributable  to  the  fluid  drag  was  the  only  mechanism  provided  by  the  interstitial  fluid. 

The  results  of  this  work  show  quantitatively  the  significance  of  size  gradations.  For  constant 
volume  concentrations,  as  the  smaller  particles  decrease  to  half  the  original  size,  the  internal  shear 
stress  could  be  reduced  by  as  much  as  ten  times.  Although  theory  checks  well  with  the  limited 
existing  data,  the  general  validity  of  the  present  theoretical  work  has  to  be  confirmed  by  future 
experimental  work  for  a  wider  variety  of  size  gradations. 

There  are  a  few  difficulties,  however,  regarding  the  description  of  particle  configurations  that 
may  introduce  errors  in  the  present  analysis.  These  difficulties  are  discussed  below. 

The  densest  volume  concentration  CQ  that  appeared  when  defining  the  average  gap  between 
particles,  and  also  in  the  final  form  of  the  stresses,  is  not  a  constant.  Nevertheless,  in  all  the  com¬ 
putations  it  is  assumed  that  CQ  is  equal  to  0.74,  which  is  the  densest  volume  concentration  for  a 
hexagonal  packing  of  uniform  spheres.  This  parameter,  CQ,  should  depend  on  the  packing  pattern 
and  the  size  ratio  of  the  two  different  particles.  When  the  small  particles  are  infinitely  small,  they 
can  fill  the  voids  between  the  large  spheres.  Hence,  in  the  limiting  case  the  densest  volume  con¬ 
centration  for  a  two-size  packing  should  be  0.74  +  (1-0.74)  x  0.74  =  0.93.  Fortunately,  values  of 
the  nondimensional  stresses  are  quite  insensitive  to  the  variation  of  CQ.  It  is  found  that  as  C  =  0.5, 
when  CQ  varies  from  0.74  to  0.92,  the  nondimensional  shear  stress  defined  in  eq  41  increases  by 
less  than  25%. 

Another  difficulty  arises  when  the  ratio  DS/DL  tends  to  be  very  small,  since  as  the  smaller  par¬ 
ticles  become  extremely  small,  a  large  number  of  small  particles  will  be  surrounded  by  small  parti¬ 
cles  only.  Namely,  as  Ds/Dt  goes  to  zero,  the  statistical  concept  used  in  estimating  the  average 
gap  between  particles  may  break  down.  It  is  believed  that  as  DS/DL  goes  to  zero,  the  average  gap 
between  small  particles  is  overestimated  by  eq  21 .  Thus,  the  actual  stresses  may  be  somewhat 
greater  than  the  theoretically  predicted  values  as  Dsi'Dl  goes  to  zero. 
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Figure  6.  Nondimensional  shear  stress  versus 
nondimensional  shear  rate;  d  =  diameter. 

(A  fter  Savage  and  Sayed  f 1 984] ,  data  from 
their  Figures  9b  and  18b.) 


Figure  7.  Particle  size  distribution  of  spheres  (after  Savage  and  Sayed  1984). 

This  prediction  is  of  the  same  order  of  magnitude  as  the  measured  data.  The  discrepancy  be¬ 
tween  the  numerical  values  of  the  data  and  the  prediction  may  be  the  result  of  the  nonuniform 
size  distribution  of  both  tests  1  and  2.  As  shown  in  Figure  7a,  the  actual  size  of  spheres  used  in 
test  1  is  not  uniform.  In  fact,  the  spectrum  is  quite  wide.  The  two-size  mixture  used  in  test  2  as 
shown  in  Figure  7b  is  not  composed  of  two  distinct  sizes  either.  However,  the  spectra  are  narrow¬ 
er  than  the  one  in  test  1 .  Based  on  the  present  theoretical  study,  it  is  expected  that  if  a  true  uni¬ 
form-size  material  is  used  in  test  1 ,  measured  stresses  would  be  greater  than  the  values  given  in 
Figure  6.  Therefore,  the  ratio  given  in  eq  42  could  be  closer  to  that  given  in  eq  43  than  it  appears. 

Since  data  points  for  only  one  size  ratio  RD  and  one  concentration  ratio  RC  are  available,  the 
general  validity  of  the  present  theoretical  work  awaits  further  experimental  study.  The  results 
of  these  experiments  will  provide  data  that  cover  the  full  range  of  the  size  and  concentration  dis¬ 
tributions. 


cf 

1 

A 


CONCLUSION 

This  report  derives  the  stresses  generated  in  a  highly  concentrated  shear  flow  of  spherical 
solids  with  two  sizes.  The  stresses  are  considered  to  be  caused  predominantly  by  binary  colli¬ 
sions  between  adjacent  particles.  The  interstitial  fluid  effect  is  neglected.  The  fluid  effect  could 
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Figure  5.  Particle  size  effect  on  shear  stress. 
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•  [2 s/2  ( RD +  RC2RD-l3li)  +  4  Ra-RD^  +  RC(fUl)l  ] 
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As  shown  in  Figure  5,  the  nondiinensional  shear  stress 


C  P%D{(~  )  (1+c)  (0.053  +  0.081  p)  " 


is  plotted  against  the  two  parameters  RD  =  DS/DL  and  Cs  for  a  typical  total  concentration  C. 


COMPARISON  WITH  EXPERIMENTAL  DATA 

Savage  and  Sayed  (1984)  experimentally  studied  the  stresses  generated  in  a  couette  device 
filled  with  dry  plastic  beads.  The  results  of  two  tests  are  shown  in  Figure  6.  Test  1  used  nearly 
uniform  beads  of  1 ,32-mm  mean  diameter.  Test  2  used  beads  of  a  7:3  mixture  of  mean  diameters 
1 .65  mm  and  0.6  mm,  with  a  weighted  mean  of  the  two  sizes  of  1 .34  mm.  The  size  distributions 
of  these  two  tests  are  shown  in  Figure  7.  The  straight-line  fit  of  the  two  sets  of  experimental 
data  shown  in  Figure  6  demonstrates  that 

t(  two  sizes)  __  I 
T<uniform)  5 

For  this  given  size  distribution,  the  theoretical  prediction  using  eq  39  gives 
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and 
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with  s  given  by  eq  21. 

If  we  use  the  results  obtained  for  v'L  and  Vg  and  substitute  eq  18,  21 , 26,  28, 29  and  32  into 
eq  2  and  3,  the  constitutive  equations  for  a  granular  mixture  of  spheres  with  two  different  diam¬ 
eters  are 
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(l+/?£»3)'^ 

in  which  /?C  =  Crs/GL  and  flZ)  =  DS/DL  ,F  is  defined  in  eq  35.  Since 


Gs  Ns  Cs  Dl, 

RGu£--nl--cl<d?  ~RC(RD) 


(37) 


(38) 


in  which  RC  =  CS/CL ,  eq  36  and  37  can  be  rewritten,  respectively,  in  the  following  more  conven¬ 
ient  forms: 


(l-t-e)3/2  (0,053  +  0.081  nfh 
.l~e2  y(l+e)  y2(l+e)2 . % 
‘4  rr  '  4  1 


[l+(/?(7)f?Z)'3]  c’/3 
2|l+(/?CV?D'4)  (C01/3  -Cl/3) 


.  F*12  .  1 

[l+(/?C)/?Z)’9/2]2  l+/?c 


(39) 


V, 

V 
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Ai7  _  nuh  2(l+e)  , 

^LL2  "  "PS  g  2  V2  KL 


AM 


”D3S 


du 


SSI  =  Ps  (1  +  €)  (0.053  +  0.081  fi)j^  (Ds+s) 


AiJ  _  nDs  2(1 +e)  /T  , 
^SS2  ~  "Ps  a  15  v2  I's  . 


(28) 


The  momentum  transfers  between  particles  of  different  sizes  are  obtained  similarly.  We  can  derive 
A/Wj_s  by  considering  a  stationary  large  particle  being  hit  by  a  small  particle  and  statistically  aver¬ 
aging  the  post-collision  momentum  of  the  large  particle.  With  the  details  shown  in  Appendix  A, 
the  results  for  A MLSl  2  and  AA/SL1  2  are 


nDl 


du  Ps+dl 


2D\ 


^sli  =  AMLS1  =  ps  (1  +  e)  (0.053  +  0.081  p)jg  (- ~^  +  s) 

*Dl  4  (l+e)D| 


AAfSL2  AA/Lg2 


"Ps  c.  _2 


s  6  7 rJ  Dl+Dl 


(29a) 

(29b) 


Average  energy  dissipation  in  a  pair  of  spheres-^LL,  E^g, 

The  average  energy  loss  in  a  pair  of  colliding  particles  with  identical  size,  as  derived  by  Shen 
and  Ackermann  (1982),  becomes 


irDl 


'LL  or  SS  P s  ^ 


•2 

V 

L  or  S 


rl-e2 


+  /i(l+e) 

7 r 


P2(  l+g)2! 
4  J 


(30) 


The  average  energy  loss  in  a  pair  of  colliding  particles  with  different  sizes  can  be  obtained  similarly. 
With  the  details  given  in  Appendix  B,  the  result  is 


CONSTITUTIVE  EQUATIONS 


The  functional  forms  that  appear  in  eq  7-1 1  are  quantified  explicitly  in  eq  18  and  27-31.  When 
we  substitute  eq  1 8  and  27-31  into  eq  6,  together  with  the  equation  for  the  equipartition  of  energy, 
eq  12,  v'L  and  v's  are 


,  2 
L 


du 


K. 


k2f 


*-*<£>■ 


L.  3 


(32) 


in  which 


.V-W 


where  A'L  is  the  number  of  large  particles  in  a  unit  volume  and  /Vs  is  that  of  small  particles.  The 
total  number  of  collisions  between  identical  particles  is 


N, 


LL 


"lK.l 

2 


(25a) 


N, 


^ s  N Is 


ss 


(25b) 


since  all  collisions  between  like  particles  are  counted  twice.  Also,  the  number  of  total  collisions 
in  the  flow  field  is 


^LL  +  ^LS  +  ^SS 


(26) 


If  we  solve  for  the  seven  unknowns-A^ ,  NLS,  Nss,  A^L,  A^s,  AgL  and  A^ -among  the  seven 
equations  (eq  22-26),  the  following  results  are  obtained: 


Ar. 


LL 


N.  — —  — 
lGVT  * 


L  ^L 
2s 


(27a) 


in  which  GK7’  =  +  GsvS’ 


r  t 

glvl  vs 

A/  =  A/  — 

^LS  "s^y  , 

(27b) 

csvs  vs 
"ss  "sGVT  2s 

(27c) 

N*  =GjAVl 

LL  GVT  s 

(27  d) 

p  _GSVS  VL 

GVT  s 

(27e) 

VS  -G^  ^ 

GVT  s 

(270 

=  GSVS  *8 

^SS  GVT  J  • 

(27g) 

Average  momentum  transferred  between  pairs  of  spheres- 
A^ll  ,  LS ,  AAfSL  >  AAfss 

When  two  particles  of  like  size  collide,  the  average  momentum  transfer,  as  determined  by  Shen 
and  Ackermann  (1982),  can  be  described  as 

lli  =  Ps ^  (l+e)  (0.053  +  0.08 \  n)  ~  (0L+s) 


8 


V"'  . 


*  i  • 


yv 


namely,  a  particle  plus  its  immediate  neighbors,  is  (VL  [(/?L  +  £>s )/2)  +  GS(DS)  for  a  cell  centered 
around  a  small  particle.  For  a  cell  centered  around  a  large  particle,  the  radius  is  GLDL  +  Gs  [(Z)l 
+  Z)s)/2] .  In  a  looser  state,  as  the  volume  concentration  becomes  C ,  the  immediate  neighboring 
particles  are  all  moved  a  distance  s  away  front  the  center  particle.  Hence,  the  corresponding  radii 
for  these  "cells”  become  GL  [((/>L  +  Ds)l 2)  +  s]  +  GS{DS  +  s)  and  6L(DL  +  s)  +  CS[((Z)L  +  Ds)l 
2)  +  s] .  The  mean  free  path  s  can  therefore  be  related  to  the  relative  concentrations: 


C_  _  +  gS^S)3  _  , 

C0  =  (GlRl  +  GsRs  +  s)3  ;  L-  L- 


+  G, 


Dl  +  Ds 


D ,  +DS 

=  “T  ~ 


+ 


(20) 


Hence 

-  Co  1/3 

s  =  /)[(-^)  -  1]  with/)  =GlRl  +  (7s/?s. 


(21) 


Frequency  of  collisions  between  spheres— , 

^LS >  ^SS >  ^LL  >  ^ls >  ^SL >  i<VSS 

Consider  a  large  particle  moving  in  a  granular  flow  as  shown  in  Figure  4.  This  large  particle  is 
being  struck  by  the  neighboring  small  and  large  particles.  iV^s  is  the  total  number  of  collisions 
per  unit  time  this  large  particle  experiences  from  the  surrounding  small  particles  and  A/^L  is  the 
total  number  of  collisions  per  unit  time  this  large  particle  experiences  from  the  surrounding  large 
particles.  The  collision  frequency  is  proportional  to  both  the  number  density  of  the  surrounding 
particles  and  their  fluctuating  speed;  hence 


irs  vs 

^LS  =ALL^lV7 


(22) 


where  Gs  and  GL  are  the  number  percentages  of  small  and  large  particles  in  the  flow  field.  Simi¬ 
larly,  for  a  fixed  small  particle  in  the  flow. 


_  *  GS  VS 

Ass=  slcl^' 


(23) 


The  total  number  of  collisions  per  unit  time  between  large  and  small  particles  in  a  unit  time  is 


^LS  ^LS 


(24a) 


or 


^IS  ^S^SL 


(24b) 
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Figure  4.  A  large  sphere  being  struck  by 
neighboring  spheres. 


vv 


hD 2  _  1  7r(Z>  sin^)2 


k™:  =  l  r 

4  ,  w 

(?)  Jo 


where  D  =  DL  oi  Ds,x2  =  D/ 2  costp  and  ^  is  defined  in  Figure  1 .  When  solved  for  K,  eq  13  shows 
that  the  average  area  of  a  sphere  cut  by  a  random  unit  surface  area  is  (n/6)D2 ,  where  D  =  Dt  or  £>s. 
The  total  areas  occupied  by  large  and  small  particles  on  a  unit  area  are  therefore  ph(nj6)D ^  and 
ps  (ir/6)Dl  respectively. 

Let  NL  be  the  number  of  large  particles  in  a  unit  volume  and  7VS  be  the  number  of  small  parti¬ 
cles.  Assume  that  in  a  unit  volume  the  number  percentage  of  large  particles  is  GL  and  that  of  the 
small  particles  is  Gs.  It  is  easily  seen  that 

nl  Ns 

«A,C‘  '  04» 

The  volume  concentration  of  solids  is 


C  =NL-~  +  . 


The  number  concentration  of  large  and  small  particles  is  solved  using  eq  14  and  15,  resulting  in 


*(GLDl+GsD3s) 


6  (GlDl  +  GsD3s) 


In  a  completely  random  situation,  the  volume  concentration  is  equal  to  the  area  concentration  of 
particle  cross  sections  occupying  a  unit  area;  hence 


nDj  7 iD\  irDl 

Cl  s»L-raPL-J-  ’  Cs  =Ns  g 


P  L  =  NLDL 


B  cgldl 
1(gldI  +  gsdd 


-  Ps  ~  »sDs 


CG  sP% 

~(GlD^GsDD 


It  is  convenient  to  now  derive  another  useful  kinematic  variable,  the  average  gap  s  between  ad¬ 
jacent  particles  as  shown  in  Figure  1 .  Let  C0  be  the  concentration  corresponding  to  a  densest 
state  for  some  fixed  proportion  of  large  and  small  particles.  In  this  state  all  adjacent  particles  are 
touching.  In  a  unit  volume  the  number  percentage  of  large  particles  is  CL  and  that  of  small  ones 
is6's.  By  definition. 


6V ,  ft.// £l  ,  Ik 

in 


t;s»— l/  —  *  — 

K\f*  W 


In  the  densest  state  a  given  particle  is  touched  by  a  number  of  large  and  small  particles;  the  per¬ 
centage  of  large  ones  is  G’L  and  that  of  the  small  ones  is  G's.  Hence,  the  average  radius  of  a  “cell," 
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